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Abstract 

Starting with any nondegenerate triangle we can use a well defined interior point of the 

^^ triangle to subdivide it into six smaller triangles. We can repeat this process with each new 

^^ triangle, and continue doing so over and over. We show that starting with any arbitrary triangle, 

-V.1 the resulting set of triangles formed by this process contains triangles arbitrarily close (up to 

^^ similarity) any given triangle when the point that we use to subdivide is the incenter. We 

^ also show that the smallest angle in a "typical" triangle after repeated subdivision for many 

*~~5 generations does not have the smallest angle going to 0. 
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1 Introduction 

Given a triangle and an interior point of the triangle we can divide the triangle up into six smaller 
triangles, also called daughters, by drawing hue segments (or Cevians) from the vertex through the 
interior point to the opposite side (see Figure [I]). The process can then be repeated with each new 
triangle with its own corresponding interior point and again repeated over and over. When the 
interior point is the centroid this corresponds to barycentric subdivision. 
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C^^ Figure 1: How to subdivide the triangle given an interior point P. 



Stakhovskii asked whether repeated barycentric subdivision for a starting triangle is dense is 
the space of all triangles, i.e., the space of triangles up to similarity where two triangles are e-close 
if the maximum difference between their corresponding angles is less than e. This question was 
answered in the affirmative [2]. 

Theorem 1 (Barany-Beardon-Carne). Successive subdivisions of a non- degenerate triangle using 
the centroid point contains triangles which approximate arbitrarily closely (up to similarity) any 
given triangle. 

Moreover, it was shown that almost all triangles became "flat" (in the sense that the largest 
angle approaches vr for most triangles as the number of subdivisions increase). Ordin \T7\ extended 
their results and showed that this also holds if we choose the interior point of the triangle to be 
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PqA + piB + P2C where A, B, C are the vertices with po + pi + P2 = 1 and pi > (the centroid 
corresponds to po = Pi = P2 = g)- 

Our main result is to show that a similar statement holds if we choose the interior point to be 
the incenter, which can be found by taking the intersection of the angle bisectors. 

Theorem 2. Succesive subdivisions of a non-degenerate triangle using the incenter point contains 
triangles which approximate arbitrarily closely (up to similarity) any given triangle. 

We will see however that there is a difference in the behavior in that not almost all triangles 
become flat as in the centroid case. 

We will proceed as follows. In Section [2] we will give a quick sketch of Theorem [T| while in 
Section [3] we will give a proof of Theorem [2] and establish several properties about this subdivision. 
Finally, in Section |4] we will give some concluding remarks. 

2 Subdividing using the centroid 

In this section we give a quick sketch of the ideas behind Theorem [T} The method of Barany et al. 
[2] was to first associate triangles with points in the hyperbolic half plane, namely each triangle T 
is associated with (up to) six points z in the hyperbolic upper half plane HI as shown in Figure l2| 
This is done by placing soTne edge of T with vertices at 2; = and z = 1 and the third vertex is 
located at the complex coordinate z with positive imaginary part. Observe that reflecting z across 
the three circles 5R(z) = ^, \z\ = 1, and |-z — 1| = 1 induces a natural action of ^3 on the hyperbolic 
half plane HI in which all six orientations of T occur. 




Figure 2: Triangles as points in the hyperbolic plane. 

Now, the centroid point of a triangle with vertices at 0, 1 and z is the point (z + l)/3, and so 
one of the corresponding daughters becomes 2(2; + l)/3 when normalized. The argument reduces 
to showing that the group of automorphisms of H generated by the map B{z) = 2{z + l)/3 and 
the above 5*3 action is dense in Aut{M), in particular for any starting z (i.e., any intitial triangle 
T) the set of all resulting points is dense in EI (i.e., dense in the space of triangles). 

Further, using results of Furstenberg j^, it follows that almost all random walks formed from 
products of B{z) and elements of 5*3 tend to inflnity (in the hyperbolic plane) as the length of the 
product increases. This then implies that almost all of the nth generation daughters have smallest 
angle tending to as n increases. By different techniques, Robert Hough [8] was able to show that 
the largest angle approaches vr and moreover was able to give asymptotic bounds for the proportion 
of triangles with angles near vr. 



3 Subdividing using the incenter 

The important step in the proof of Theorem [T] was to find a way to associate triangles with points 
where the action of finding a daughter triangle was natural. The first step in proving Theorem [2] is 
to do the same. However, we will find it more convenient to associate each triangle with a point (s) 
in M^ where the coordinates are the angles. The set of all possible triangles (including degenerate 
cases), denoted P, is the intersection of the hyperplane x + y + z = tt with the first octant (see 
Figure pi) . Note that P also corresponds to a two dimensional equilateral triangle (see [H [5l [TOl [TT] 
for previous applications involving P). 
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Figure 3: Representing triangles as points on P in M^. 

As noted in the introduction the incenter is found by the intersection of the angle bisectors. So 
the angles of the new triangles created by subdivision are linear combinations of the angles of the 
original triangles (hence the reason it is more convenient to work with P). In particular, if we let 
t = (a, /?, 7)* denote a triangle then the six new triangles are found by Mjt where 
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Figure 4: The image of P under the six maps Mj. 



Observation 1. The union of the MiP covers P. 

This can be seen by examining Figure |4j Alternatively, this says that every point in P has a 
preimage in P under some Mi. For the case that t = (a,/3,7)* with a < /? < 7 we have that 
(Mi)^^t = (2a, 2/3 — 2a, a — /3 + 7)* is a preimage of t in P. Other possible arrangements for the 
ordering of a, /3, 7 can be handled with the remaining Mj. 

Observation 2. If we let \\-\\ denote Euclidean distance, then ||Mj(t — s)|| < 
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To see this we can put P into M^ by t = (a,/3,7)* — )• ( /^ ; q) j note that this will preserve 

distance. Under this map we would also have that Mit — )• ( ° J^ , ^) . If we now let s = (a', /3', 7') 
then a calculation shows 

-lit - slP - llMit - Misf = -(/?- P'f > 0. 

The result now follows for Mi and similar calculations establish it for the remaining Mj. 

We now prove Theorem [2j Let q be the initial triangle we apply subdivision to. We need to 
show that for any triangle t and e > there is some sequence of ij so that 

||(Mi,Mi2---Mi,q)-t|| <e. 

Choose k sufficiently large so that 7r(^) < e. By Observation 1 we can successively find a /cth 
generation preimage of t in P, which corresponds to multiplying by an appropriate (Mj)~^ at each 
step. Denote this preimage by 

(M,J-i-..(M,,)-i(M,J-H, 

(where the ij are chosen according to how we construct the preimage). Repeatedly using Observa- 
tion [2] we have 

||(Mi,Mi,---Mi,q)-t|| = ||Mi,(M,2---Mi,q-(MiJ-it)|| 

< ^||iV4---M,,q-(M,J-H|| 

< ... 

< <^) < e- 

In the last step we used that points in P are at most distance vr apart. This finishes the proof of 
Theorem [2l 

The limiting distribution 

In fact more can be said about the iterated subdivision of triangles using the incenter. Namely, 
since the maps Mj are contracting with Lipschitz constant v3/2 then it follows (see [6j) that there is 
a fixed limiting distribution on P that the process converges to. Further it converges exponentially. 
To get some sense of what this limiting distribution looks like we can simply start with any 
triangle (in our case we will use an equilateral triangle) and plot all of the nth generation daughters 
for some n in P. This is done for n = 5 in Figure [5a| 



Examining Figure 5a we see that the daughters seem to fill in most of P (agreeing with Theo- 
rem 2|. However, a patient count will reveal that there are far fewer than 6^ triangles in Figure 5a 



This is because some points have been mapped onto several times (a consequence of starting with 
such a symmetric triangle) . So to get a better sense of the limiting distribution instead of plotting 
the individual triangles in P it is better to look at a histogram. We will divide P into a large 
number of small regions and then shade each region according to the number of triangles that fall 
into that region, the darker a region is the more triangles fall into that region. In Figure 5b we 
give the histogram for n = 12 generations starting with the equilateral triangle. 
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Figure 5: The distribution of the nth generation daughters using the incenter. 



Very httle is known about the hmiting distribution. ExperimentaUy, it appears that the densest 



point on the hmiting distribution (i.e., the darkest region in Figure 5b) corresponds to the triangle 



(|, |, ^)*. This is hkely because this triangle corresponds to an eigenvector of eigenvalue 1 of 
two of the Mi. In other words under subdivision using the incenter this triangle has two of its 
daughters which are similar to it (see Figure ph. No other triangle has this property, and the 
triangle (if , if , if )* is the only other triangle with one of its daughters similar to itself, but this 
triangle does not appear to play a significant role in the distribution. 




Figure 6: The triangle 
to the original). 
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* subdivided using the incenter (the shaded triangles are similar 



If we were to draw the histogram for n = 20 or n = 50 and compare it to Figure 5b we would see 
almost no perceptible difference between them. This is because, as we noted above, the convergence 
to the limiting distribution is exponential. Or put another way, if we look at what happens when 
we map a triangle under n applications of the Mi in P then knowing what the last few steps that 
we applied gives us a good handle on where we are in P. This is essentially the heart of the proof 
of Theorem [2j 



For example, the region M1M2M3P corresponds to a triangle with vertices in P at (f 
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(|, J, ^)* and (|, |, ^)*. In particular, looking at all the daughters for n large at least 1/6^ of 
the daughters must lie in this subregion of P (i.e., 1/6^ of the possible products of the Mj will 
have M1M2M3 as the leading term). Since points inside this subregion of P must have minimum 
angle at least vr/S then we have that at least 1/6^ of the daughters in the nth generation must have 
minimum angle at least vr/S. 

Of course, by looking at larger products and looking over more prefixes we can say a lot more 
about what happens with the minimum angle. For example, in Figure [7a| |7b| and [7c] we have 
plotted all of the triangles that are the boundary of P under all possible second (MjMjP), third 



(MjMjMfcP) and fourth {MiMjlSd^M^Y*) maps of P into itself. (The number of such triangles is 
large so it is difficult to pick out the individual triangular regions. It also is interesting to note the 



similarity with the histogram in Figure 5b 






(a) Second generation (b) Third generation (c) Fourth generation. 

Figure 7: Resulting images of P for different generations. 

We can now bound the limiting cumulative distribution function (CDF) for the smallest angle 
in the limiting distribution of triangles. This is done by considering all resulting 6" images of P in 
the nth generation. Then for any angle 6, a lower bound for the number of triangles with minimum 
degree 9 (or less) is found by counting the number of the images of P which have largest minimum 
degree at most 9. Similarly, an upper bound for the number of triangles with minimum degree at 
most 9 (or less) is found by counting the number of the images of P which contain a triangle with 
minimum degree at most 9. Doing this for n = 11 gives Figure pj (By comparison the limiting 
CDF under subdivision using centroids is the constant function 1, showing that these two methods 
of subdividing are fundamentally different.) 
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Figure 8: Upper and lower bounds for the CDF for the smallest angle in the limiting distribution. 



4 Conclusion 

We have seen that like the centroid, when doing subdivision using the incenter the resulting triangles 
are dense in the space of all triangles. However, unlike the centroid the smallest angle in a typical 
triangle does not tend to 0. This is important since certain methods can fail when the subdivision 
creates a large number of triangles with minimal angles going to as n gets large (see O [THl [T^ ) . 
One interesting question is to understand the limiting distribution for repeated subdivision 



using the incenter, an approximation of which is shown in Figure 5b 



One can also consider what happens for subdivision using other interior points. For example, 
the Gergonne point is found by taking the inscribed circle in the triangle and connecting a vertex 
to the point of tangency on the opposite edge; these three lines intersect at the Gergonne point. 
When using the Gergonne point to subdivide it is known [Ij that the triangles are not dense in 
the space of all triangles. In Figure [9a] we have given a histogram of P for the tenth generation of 
subdividing using the Gergonne point (notice the large white spaces where there are no triangles). 

An interesting point for which little is known about what happens after repeated subdivision 
is the Lemoine point, which is found by taking the lines from a vertex through the median of 
the opposite edge and then flipping them across the angle bisectors; these three lines intersect at 
the Lemoine point. In Figure 9b we have given a histogram of P for the eleventh generation of 



subdividing using the Lemoine point. It is currently unknown if this method of subdivision is dense 
in the space of all triangles and what the limiting behavior is (there is some experimental evidence 
that the triangles become flat, but the convergence seems to be relatively slow). 

More information about the Gergonne and Lemoine points, as well as a large number of other 
interesting points available to investigate, can be found online (see [9]). More information about 
what happens under repeated subdivision using a central point can be found in [5]. 
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(a) Gergonne point (b) Lemoine point 

Figure 9: Histograms for the Gergonne and Lemoine points. 
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